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Abstract. We present the exact solution of the Richardson-Gaudin models 
associated with the SU(3) Lie algebra. The derivation is based on a Gaudin 
algebra valid for any simple Lie algebra in the rational, trigonometric and 
hyperbolic cases. For the rational case additional cubic integrals of motion are 
obtained, whose number is added to that of the quadratic ones to match, as 
required from the integrability condition, the number of quantum degrees of 
freedom of the model. We discuss different SU(3) physical representations and 
elucidate the meaning of the parameters entering in the formalism. By considering 
a bosonic mapping limit of one of the SU(3) copies, we derive new integrable 
models for three level systems interacting with two bosons. These models include 
a generalized Tavis-Cummings model for three level atoms interacting with two 
modes of the quantized electric field. 
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1. Introduction 

The Richardson-Gaudin (RG) exactly solvable models [HE] can be traced back to the 
exact solution of the BCS Hamiltonian given by Richardson in the early 1960's [3] 
and to the integrable spin model developed by Gaudin in the seventies [I]. They are 
diagonalizable by Bethe Ansatz techniques and are the so-called classical limit of two 
dimensional vertex models. We refer the reader to [5], where the connection between 
the RG models and the inhomogeneous XXX vertex model with twisted boundary 
conditions is established in detail for the case of the rank-1 SU(2) algebra. Previously, 
in [B] it was shown that the solution of the Gaudin models [those without linear term 
in the integrals of motion, cf. (|19[) below] associated with more general Lie algebras, 
can be used to get a solution of the corresponding Classical Yang-Baxter Equations 
(CYBE). This connection has been largely exploited to diagonalize the RG models 
(see EJ [7] for some examples) . In practice this method consists of using the already 
known solutions of the general-Lie-algebra Yang-Baxter Equations (obtained by using 
the Quantum Inverse Scattering method) to obtain the corresponding solution of the 
RG models by taking the classical limit of the respective Bethe equations. From a 
Conformal Field Theory context, Asorey, Falceto and Sierra [8] obtained the integrals 
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of motion and respective eigenvalues of the rational RG models for any simple Lie 
algebra as a limiting case of the Chern-Simons theory. 

Here we follow an alternative and more direct approach to diagonalize the RG 
models. Though essentially equivalent, it differs in practice. This approach, where 
no reference to the YBE is necessary, is similar to the presented by Ushveridze [9] for 
the rational case. The method is based on the introduction of an infinite dimensional 
algebra (the Gaudin algebra) associated with the Lie-algebra of a simple group. By 
taking a Casimir-like operator in the Gaudin algebra, one gets the transfer matrix 
of the YBE approach and, consequently, a set of independent quadratic integrals of 
motion. In [lOj this formalism was extended to the trigonometric and hyperbolic cases 
for the rank-1 algebras SU(2) and SU(1,1). In this contribution we extend the Gaudin 
algebra to the trigonometric and hyperbolic cases for any simple Lie algebra. 

As an application of the method, we work out in detail the RG models associated 
with the rank-2 SU(3) algebra. Most of the physical applications of the RG models 
presented so far are based on the rank-1 algebras SU(2) or SU(1,1). They cover a wide 
variety of physical problems ranging from pairing Hamiltonains for fermion or boson 
systems [HIE], to spin models or generalized Tavis-Cumming models [IB]- Recently, 
some applications to physical problems of RG models based on higher rank algebras 
have been published. The detailed derivation of the exact solution was considered for 
the SO (5) algebra of proton-neutron pairing in |14j and with more generality including 
numerical applications in 15J. In [16] the RG model based on the non-compact 
SO(3,2) algebra is used in the context of the Interacting Boson Model II, which 
describes the interaction between two different species of bosons (it and v) for each 
non-degenerated level. In [7] the SU(4) RG models are used to describe the interaction 
between two superconducting systems. Here we discuss different physical scenarios of 
potential interest, associated with the SU(3) RG models. They include dipole-dipole 
interactions between three level atoms, isospinorial pairing and generalized Tavis- 
Cummings models of three level atoms interacting with two different bosonic modes 

Contrary to rank-1 RG models, in higher rank algebras the number of independent 
integrals of motion is not exhausted by those obtained from the quadratic Casimir- 
like operators. More integrals of motion can be obtained by considering higher degree 
operators. Unfortunately, as it was shown in |18j . the Casimir-like analogy to obtain 
integrals of motions is valid up to degree-three Casimir operators. More general 
formulae are needed to get integrals of motion of degree greater than three |19j . 
Nevertheless, since the SU(3) algebra has two independent Casimir operators of degree 
two and three respectively, the Casimir-like analogy can be used to obtain the complet 
set of integrals of motion. In this contribution it is verified that, for the rational SU(3) 
RG models, the number of independent integrals of motion coincides with the number 
of quantum degrees of freedom as defined in [20]. For other higher-rank algebras, 
which have at least quartic Casimir operators, the more general formulae of [19] are 
needed to obtain the complete set of integrals of motion. 

This paper is organized as follows. In section two, the general Gaudin algebra 
and quadratric Casimir-like operators are introduced. A particular realization of the 
Gaudin algebra in terms of a direct product of L Lie algebra copies and complex 
valuated functions is introduced. The conditions to be satisfied by these functions are 
established, which are a generalization of the Gaudin conditions of the SU(2) case. 
Three particular solutions to the Gaudin conditions are found (rational, trigonometric 
and hyperbolic) and it is shown how to obtain the corresponding sets of quadratic 
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integrals of motion from the Casimir-like Gaudin operator. Special interest is paid to 
the introduction of a linear term in the integrals of motion through the addition of 
a constant shift to the Cartan members of the Gaudin algebra. Explicit formulae for 
the quadratic integrals of motion and their eigenvalues are given. In section three, 
focusing on the SU(3) algebra, we apply the formalism of section two to obtain closed 
expressions for the quadratic integrals of motion and their eigenvalues in the more 
general scenario of arbitrary SU(3) irreducible representations. As discussed above, 
in order to satisfy the condition of integrability more integrals of motion are needed, 
which must be at least cubic in the generators. For the rational version we present new 
sets of integrals of motion coming from the Casimir-like Gaudin operator of degree 
three, and it is verified that the total number of quantum constants of motion coincides 
with the number of quantum degrees of freedom. In section four, the physical meaning 
of the variables appearing in the integrals of motion and eigenvalues is discussed for 
different SU(3) physical realizations. Additionally, the limit of infinite degeneracy for 
a copy of the SU(3) algebra is used to get a family of integrals of motion related to a 
generalized Tavis-Cummings model for three level atoms interacting with two species 
of bosonic excitations. Conclusion are given in the last section. 



2. Richardson- Gaudin models from a Generalized Gaudin algebra 

Let us begin by considering a simple Lie algebra, expressed in terms of its Cartan- Weyl 
decomposition [2T1I22] : 

[S a ,S fc ] =^C c afc S c , (1) 

c 

where the non-zero structure constants are given by: 

9 

a - ~ U a - a ! °i - |^|2 a ' °"+/3 ~~ Q >' 3 ' 

the Latin index runs over the r Cartan- Weyl members of the Cartan subalgebra 
(r =rank of the group), and Greek indexes refer to the roots, a 1 denotes the component 
i of the root a. 

Associated with this algebra we propose the following infinite dimensional Gaudin 
algebra: 

[S a (A), S 6 ^)] = J2 Cf{X b (X - M )S C (A) - X* a (\ M )S c (m)}, (2) 

c 

where X a (A) are meromorphic functions associated with the member a of the Gaudin 
algebra and A is a complex variable. From these commutation rules, it can be proved 
that the Casimir-like operators of the Gaudin algebra, 

K(A)^ 5ah S a (A)S b (A), (3) 

ab 

where g a b is the inverse of the Killing form, commute among themselves: 

[K(A),K(/x)]=0 V A, [i. (4) 
A realization of the Gaudin algebra is given by: 

L 

S a (A)= £S£,X (z m -A), (5) 
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where the index m runs over L different copies of the Lie algebra, SJ^ is a generator of 
the m-th copy, and z m is a set of completely free real parameters. The operators S a (A) 
act upon the space Vj. <8> ... ® Vl, with V^j an irreducible representation (irrep) of the 
Lie algebra. Given this specific realization, the commutation rules @ impose some 
conditions on the functions X a (X). All the functions associated with the elements of 
the Cartan subalgebra are equal, real valuated and antisymmetric: 

Xi(X) = Z(X) = Z*{\) Vt (6) 
Z(X) = - Z{-\). (7) 

The reality condition on the functions Z(X) comes from keeping the hermiticity of the 
Cartan-subalgebra in its Gaudin counterpart [S J (A) = (S l (A))^]. On the other hand, 
since S Q (A) = (S~ Q (A))t the functions associated with the Gaudin ladder operators 
must satisfy: 

X a (X) = X*_ a (X) Wa. (8) 

Additionally they have to be anti-hermitic: 

X a (X) = -X*(-\). (9) 

The last conditions that have to be satisfied by the functions Z{\) and X a (A) are a 
generalized version of the S[/(2)-Gaudin conditions: 

X a (n - X)X a (X) + Z{\ - n)X a (p) - Z{\)X a {n) =0 V a 

Xp(n-X)X a+f3 (X)+X a (X-fi)X a+p ( i i)-X a (X)X p (ji) = Va,/3 such that 

a + (3 is a root. (10) 

Using the structure constants, eqs. @, ^ and ([5]), and denoting the Gaudin-operators 
by the more common notation S*(A) = H'(A) and S Q (A) = E Q (A), the following form 
for the Gaudin algebra is obtained: 



[ff(A),E«( M )] =a l {X a (X- f i)E a (X)-Z(X-n)E a (n) } (11) 

[E«(A), E-«( M )] = j^X* a (X - fi){a ■ H(A) - a ■ H( M )} (12) 
\ a \ 

[E«(A),E^(/i)] =N a>l3 {X (X-riE<*+0(X)-X*(X-riV°'+t 3 (ri}, (13) 

where a ■ H = ^[ a*H\ From equation (fT2"]l it is clear that the Cartan subalgebra 
members of the Gaudin algebra are defined up to a constant term, which can be freely 
added without altering the commutation relationships: 

ff(A)=^Hj 7l Z(z m -A) + fc 4 l, (14) 

m 

where 1 is the unity operator in V\ ® ... ® Vj,, and k l are r free real parameters. 

In the rest of the paper, we will only consider the case where all the functions 
associated with the positive root operators are equal (X a (A) = X(X), Va > 0) and 
let for the future a more general discussion. In this case, three solutions to the Gaudin 
conditions (fTU)) are: 

• The rational solution 

X(z n — z m ) = Z(z n — z m ) = (15) 
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• The trigonometric solution 

\ = i 

sm(z„ - z 



exp[i(z„ — z m )\ 

X(z n - z m ) = — — — , Z(z n - z m ) = cot [z n - z m )(16) 




• The hyperbolic solution 

X{z n -z m ) = ° xp ^ Zn ?HL- i Z(z n ~z m ) = coth (z n - z m ).(17) 

smh(z„ - z m ) 

The Casimir-like operator defined in , acts as a generator of the quadratic RG 
integrals of motion. From (fTU)) and their solutions (|15I16I17[) . it is easy to show that 
the operator K(A) can be written as: 

L „(2) L 

m—l m—1 

where C m is the Casimir of degree two of the m-th Lie algebra copy. The operators 
R m read: 

L ( I 12 

- z m )E^E-« + X*(z n - 2 m )E m Q E«) 



+Z(z n -z m )J2K l K), (19) 

i=l ) 

where the linear term £ TO = Yli=i ^H[„ is a member of the m-th Cartan subalgebra 
coming from the constant term shift introduced in (TH]) . Note that the coefficients k % 
do not have an index to, i.e., even if each £ m belongs to a different copy of the Lie 
algebra, the same linear combination respect to the corresponding basis of the Cartan 
subalgebra is assumed for all the copies. 

The commutativity of the Casimir-like Gaudin operators Q, implies 
commutativity among the set of operators R m : 

[R m ,R n ]=0 Vn,m=l,...L. (20) 

The eigenfunctions and eigenvalues of the set of integrals of motion R m are obtained by 
considering a Bcthc ansatz, which is solely written in terms of the lowering operators 
of the Gaudin algebra and a set of parameters to be determined. This calculation is 
performed in [S] for the rational case. Repeating the same steps for the trigonometric 
and hyperbolic cases is a straightforward but laborious task. It consists of applying 
the Casimir-like Gaudin operator to the ansatz. This application yields one term 
proportional to the original ansatz and terms which are not proportional. From the 
former one, one can obtain the eigenvalues of R m by performing a similar expansion 
as in (|18[) . By imposing the annulation of all the non proportional terms, one gets 
the equations that determine the parameters entering in the ansatz. We present the 
results. The eigenvalues of the operators R m are given by: 

r L r M a , , 2 

ab n (n^m) o=l fc=l 

(21) 
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where a a are the simple roots of the algebra (there are as many simple roots as the 
rank of the algebra). The r coefficients £ a are related to the parameters k l of the 
linear term through: 

e,„ = ]>> i Hj„ = ]rr^ 6 iC, (22) 

i ab 

where are the members of the Chevalley basis of the m-th Cartan subalgebra. The 
Chevalley basis is defined as 

2 

\a a \ 2 

where a a is one of the r simple roots and H m is the already introduced Cartan- Weyl 

basis of the Cartan subalgebra. F ao is a r x r matrix called the quadratic form of 

1 1 2 

the algebra, which is related to the Cartan matrix through F a b = (A^ 1 )^ . The 
Cartan matrix codifies completely the structure of the algebra and is defined in terms 
of the simple roots: A ao = , p (a a , «&), where the scalar product of the roots is 

defined as: (a a ,ctb) = J2i a a a l- are the Dynkin labels (the eigenvalues of the 
Chevalley basis) in the highest weight state of the m-th irreducible representation V m , 
and the product A m • F ■ A„ is equal to ^ Qb A^F^A*. The variables E% entering in 
(|21[) determine the common eigenf unctions of the operators R m , and are the solutions 
of the following Richardson-Bethe equations: 

r M b L 

£ E A baZ{E\, - E a k ) - X: KnZ{z m - El) - T, (k = 1, Ma)- (23) 

b=l fe'=l m=l 

The number of these parameters E% is: 

M « = rpE f 4E A ™-4 (24) 

1 a| 6=1 V m J 

where Aq are the eigenvalues of the overall operators (the sum of the Chevalley 

basis over all the copies). These overall operators commute with the integrals of motion 
R m (see ([33)) below), therefore their eigenvalues (Aq) are conserved quantities for any 
Hamiltonian defined as a function of the integrals of motion R m . Equations (fT9)) . (|2~T|) . 
and |23|) extend the results presented in [HE] for the rational case, to the trigonometric 
and hyperbolic ones. 



3. The SU(3) algebra 

3.1. Quadratic integrals of motion 

In this section we will apply the formulae of the previous one in the particular case of 
the rank-2 algebra SU(3). We consider L copies of a C/(3) algebra: 

[K a/ g TO ,K a ' / g/ TO ] = Sp a iK. a fjt m — 5 a pi~K. a i p m , (25) 

where the index m refers to the m-th copy of the Lie algebra, and a, /3 = 1, 2, 3. From 
these commutations rules it is easy to prove that [Kn m + K22m + K33 m ,K ai a m ] 
Va, (3. Therefore this sum must be proportional to the operator unity in the SU(3)- 
irrep V m and, consequently, an integral of motion: 

Kum + K.22m + K33 (26) 



SU(3) Richardson- Gaudin models 7 

This condition reduces the number of independent generators from 9 to 8, the 
dimension of the SU (3)-algebra. A Cartan decomposition of the SU(3) algebra is: 

• A maximal abelian subalgebra of hermitian operators (the Cartan subalgebra Ti.) 
is provided by the operators K. aam . However, due to the condition ([26]) only two 
of them are independent. Two different Cartan subalgebra bases adequate for our 
purposes are the Cartan- Weyl basis (we are using the normalization convention 

of my- 

ttI Kn m — K33 m 2 2K22m — K33 m — Kn m . . 

H " = 71 H " 1 " VE ' (27) 

and the Chevalley basis: 

hi = Knm — K22m n m = K.22m ~ K33 m . (28) 

These bases generate the Cartan subalgebra H m — span(H,^, H^ n )=span(h^, h^J. 

• The positive root vector space is spanned by the raising operators: 

= span(K Q/3m ) with (3 > a. (29) 

• The lowering operators are the hermitian conjugated of the previous ones. The 
negative root vector space is: 

£ ~ = span(K Q/ 3 m ) with (3 < a. (30) 

The m-th Lie algebra of SU(3) is given by the direct sum Csu(3)m = W m £+. 

The roots of the algebra can be obtained from the commutation relations between the 
Cartan- Weyl basis and the root vectors: 

[H^, K a/ g TO ] = a* a pK a p m a^f]. (31) 

The algebra SU(3) has three positive roots(/3 > a), they are: 

"12 = ("12, "12) = (1,-^3) 
"23 = ("23, "23) = T/ft 1 '^) 



"13 = {<x 13 ,a 13 ) 



(V2,0). (32) 

It is easy to see that the simple roots are ai2 and 0123 > from which all the other roots 
can be obtained by a linear combination of integer coefficients. The non-simple positive 
root is ai3 = a±2 + "23, whereas the negative roots are: 021 = —a±2, a 3 2 — —023, 
and «3i = — a 12 — a 2 3. All the roots have square norm equal to 2: |a a ^| 2 = 2. The 
/ 2 -l\ 

Cartan matrix is: A = I . From this expression we obtain the quadratic 



x _1 2 

1 ( 2 1 



form: F — ., , 

J \ 1 z 

By using all these results, we can write the integrals of motion (fT9|) for the SU(3) 
case. The rational integrals of motion read: 

K m = -eK22 m -(e+e)K 33m +c m i+ £ ^e.k^k^ (33) 
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For the trigonometric and hyperbolic cases [ct(z m — z n ) = cot(z m — z n ) and 
ct(z m — z n ) = coth(z m — z n ) respectively], the integrals of motion are: 

R m = - C 1 K 22m - (e + e 2 )K 33m + C m l (34) 

L 

ct(z n z m 



+■ E 



a 



/3>a 



The constant term appearing in (|33j) and (|34f is: C m = — (n m /3) J2i(i^ m ) Z( z i — z m )ni + 



(2£i±£), 



3. 2. More integrals of motion 

By direct calculation, it can be shown that the overall operators 

L 

N Q = 22 K aam , witha = 1,2,3, (35) 

m— 1 

commute (in the trigonometric, hyperbolic, and rational cases) with the integrals 
R m . One of these operators must be independent of the already introduced 
integrals of motion. From (|2l)|) it is clear that Ni + N 2 + N3 = Ylm n ™>> ano - 
from the integrals of motion ([33]) or ([53]) it is derived that ^2 m (C m l — R m ) = 
^ 1 N2 + (C 1 + £ 2 )N 3 . This fact implies that the integrals R m do not exhaust all 
the possible integrals of motion. In [18] it was shown for the rational case that the 
Casimir-like Gaudin operators of degree three commute among them and with the 
quadratic Casimir-like operators ©. The SU(3) Casimir operator of degree three 
is: C^ 3 ) = ^ a/37 K Q/ 3K / 3 7 K 7Q , from here it is clear that the Gaudin counterpart 

must read: K^ 3 ^(A) = ^ Q/37 K Q ^(A)K / g 7 (A)K 7Q (A), with K Q/ g(A) a generator of the 
SU(3)-Gaudin algebra. This family of operators satisfies, at least for the rational case, 
[K( 3 )(A),K( 3 )(^)] = and [K^ (A), K(fi)] = 0, V A, fj,. By expanding K( 3 )(A), as we 
did for K(A) in (JTHJl, we get: 

K (3) (A)^-V M Cm +3V / 3m -3Vi^ + VU (36) 

~ L (A - z m y ^(a-^) 2 ti x ~ z ™ h 

(3) 

where Cm is the cubic-Casimir of the m-th Lie algebra copy, and the operators p3 m 
and R,3 m are: 



3m 



L L js- TV" 



E 



, , . (z?n — Z n ) 

h ^ y A a A a i\. aam (37) 



tl-3m, = r / J A a iS~ a /3 m ]&.pa m . [6S) 
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The coefficients A a are related to the parameters of the linear term in R m : A\ = 
(2^+e 2 )/3, A 2 = (C 2 -e 1 )/3 and A 3 = -(£ x +2£ 2 )/3. The operators P 3m , R 3m and 
R m form a complete set of mutually commuting operators. From the expression for 
the integrals P3 m , it is easy to show that the three operators N Q ([55)1 can be expressed 
as a linear combination of the operators R m , n m and P3 m . The commutativity of the 
cubic Casimir-like Gaudin operators has been proved in [18j for the rational case, we 
think this result can be extended to the trigonometric and hyperbolic ones. 



3.3. Highest weight states and Quantum Dynamical Degrees of freedom 

The highest weight states (unique for any finite irrep of a simple algebra) are defined 
by the conditions: 

K Q ^ m |A m ) = for all p > a. (39) 

The eigenvalues of the members of the Chevalley basis in these states are integer 
positive numbers and allow us to label the irrep, 

h^|A m ) = A^|A m ) = (ki m - k 2m )\A m ) 

h 2 m \A m ) = A 2 m \A m ) = (k 2m - fc 3m )|A m ), (40) 

where k am are the eigenvalues of the operators K aQm in the highest weight state, and 
we have used (|28jl. The numbers [k\ m k2m k?\ determine the Young diagram of the 
irrep V m and satisfy: ki m + k 2m + fc 3m — n m . All the members of the irrep can be 
obtained by acting the lowering operators upon the previous highest weight state. 

In general, once we have completely established for each m the values of the 
set (A^,A^,n m ) (or equivalently [k\ m k 2m ^m]), we need three extra numbers for 
each copy of the Lie algebra to determine a complete basis of the quantum system. 
Consequently a complete basis for the Hilbert space V\ ® ... Vl is: 

ni ... n m ... n L \ 

|¥>= (A*, A?) ... (Al,Al) ... (Ai,A|) \ (41) 

^11,^21,^31 ■•• tlm,t 2r rnHm ■■■ ^XL^ILi^L I 

where p is a multiplicity number to take account of any other completely degenerated 
quantum number not considered explicitly here. The labels t\ m , t 2m , t% m are the three 
necessary quantum numbers to characterize completely a state within a SU(3) irrep. 
The number of non completely degenerated quantum numbers necessary to determine 
unambiguously a basis' member is the number of quantum dynamical degrees of 
freedom of the system [2CL From (|4"Tj) it is deduced that, in the present case, the 
number of quantum degrees of freedom is: 

d = 3L. (42) 

To guarantee the integrability of the system we need 3L integrals of motion. L of them 
are provided by the operators R TO , the other 2L integrals of motion are the polynomial 
operators of degree three (R 3m and P 3m ) introduced in the previous subsection. 



S.^. Eigenvalues 

We can explicitly write the eigenvalues of the operators (|33|) and (|34| , which are the 
SU(3) version of the general formula (f2Tj) : 

r m — ~ £, k 2m — (£ + £ )k^ m + C m + ^ I Z{z n — z m ) ^ k am k an J 

n(n^m) \ a=l / 
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Mi M 2 

+ (hm ~ k 2m ) ^2 Z(z m - E k ) + (k 2m - k 3m ) ^2 z ( z m - wj), (43) 
k=l 1=1 
where C m is the same constant found in the operators ([55]) and (|34[) . and we have 
redefined the parameters E k = and a->; = -E ; 2 . These parameters determine the 
eigenfunction of the operators R m , and are the solutions of the following Richardson- 
Bethe equations: 

Mi M 2 L 

2Z{E k ,-E k )-^2z(oj l -E k )-J2( k im-k 2m )Z(z m -E k )=^ (44) 

k'(k'^k) 1=1 m=l 

(fc = l,..,Mi) 

Mi M 2 L 

-^Z(Ei-ui)+ 2Z(uu-ui)- ^(k 2m -k3m)Z(z m -cj l ) =e (45) 

fc=l V{l'^l) m=l 

(l = l,..,M 2 ). 

To determine the number of parameters E k and (Mi and M 2 respectively), note 
that the members of the overall Chevalley basis can be expressed in terms of the 
integrals (f3"S"|) , 

^h.^=N 1 -N 2 , £X = N a -N 3 , 

m m 

consequently their eigenvalues are: (AJ,Aq) = (JVj. — N 2 ,N 2 — N 3 ). By using this 
result, the relation Yl m {kim + k 2m + k 3m ) — Ni + N 2 + N 3 , the quadratic form, and 
the labels of the SU(3) irreps (f^D|) in the general formula (pH)) . one gets: 

L L 

M x = kim - JVi M 2 = N 3 - ^ fc 3 ™ . (46) 

rn rn 

The values Mi and M 2 allow us to label the different invariant subspaces of the 
Hilbert space. We denote these subspaces by V{Mi,M 2 ) C Vi <8> Vr,. The different 
solutions of the Richardson-Bethe equations define a set of eigenfunctions which span 
entirely the subspace V{M\, M 2 ). By considering all the possible values of Mi and M 2 
and the set of complete solutions of the corresponding Richardson-Bethe equations, 
we get a basis for the entire Hilbert space V\ ® ... ® Vl = ® Ml M V(Mi, M 2 ). 

For a given set of L SU(3) irreps ( [ki m k 2m k 3m ] with m = 1, L), the possible 
values of Mi and M 2 are: 

< Mi < A13 
max (0, Mi - A 12 ) < M 2 < A 13 + min (0, Mi - A 12 ) , 
with A 13 = J2 m (kim ~ k 3m ) and A 12 = J2 m (kim - k 2m ). 

4. Physical models related to the SU(3) RG models 
4-1- Particle-hole representation 

An explicit realization of the U(3) algebra (|25|) can be obtained by considering Cl m 
three-level atoms of type m: 

K-af3m = ^ ] C a/i)n C /3jUm) (47) 
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Figure 1. (a) Dipolc-dipole transitions between atoms of type m and n. (b) 
Pictorial representation of a highest weight state for atoms of type m. 



where c' aflm and cp^ m are fermion creation and annihilation operators respectively. 
The indexes a and f3 label each of the three levels, whereas the index n runs over all 
the atoms of type m. We can introduce different type of atoms (let us say L different 
types), each type of atom associated with a copy of the SU(3) algebra (m — 1, .... L). 

The operators K a/ 3 m have a simple physical interpretation. For a=(3 the operators 
are the number operators: K QQm = n am . The raising operators (a < 0) take a particle 
(excitation) from a level to a lower one, whereas the lowering operators (a > (3) take 
a particle from a level to a higher one. In this context, the condition ([26]) is simply 
the conservation of the number of particles in each three level atom. 

The linear terms appearing in the integrals of motion (|33[) and (|34[) are related to 
the energies of the atoms' levels, whereas the quadratic terms represent two kinds of 
interactions: (a) the terms with a ^ (3 include three different dipole-dipole interactions 
among the atoms, which are associated, respectively, with the transitions 1 <— > 2, 1 <— > 3 
and 2^3 (see figure 1), (b) the terms with a=(3 are monopole interactions (n Qm n Q „). 

We discuss now the meaning of the labels that characterize the SU(3) irreps (the 
set (A^, A^J and n m , or [ki m , &2 m , k 3m }). Without loss of generality, we can consider 
ki m = f2 m (rt m — ki m is the number of atoms without any particle in their levels, which 
are completely decoupled of the rest and do not contribute to the r m eigenvalues). The 
atoms of each type can be classified according to the number of particle in their levels: 
(a) unblocked atoms are those with only a particle in their levels. These atoms have 
all the dipole transition allowed, (b) semi-blocked atoms are those with two particles 
in their levels, the dipole transitions between the occupied levels are forbidden by 
Pauli blocking, (c) blocked atoms are those with all their levels occupied, these atoms 
have all the dipole transitions forbidden, and interact with the others only by the 
monopole interaction. The Dynkin labels of the irreps determine the number of atoms 
of each type in the previous situations: A^ = ki m — < indicates the number 
of type-m unblocked atoms, A^ = — k^ m is the number of semi-blocked atoms 
of type m, and the value k^ m indicates the number of blocked atoms of type m. A 
representation of a highest weight state is depicted in figure 1. The case where all the 
atoms are in a unblocked situation corresponds to A^ = O m , Aj^ = and k^ rn = 0. 
Note that integrals of motion (|35[) guarantee that the total populations (irrespective 
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on the type of atom) in each level (N a , a = 1,2,3) are conserved quantities. The 
number of parameters and in the Richardson-Bethe equations Prj| has a simple 
meaning: the number of E^s (Mi) is equal to the number of atoms (of any type) with 
the first level unoccupied, whereas the number of w/'s (M2) is the overall number of 
unblocked and semi-blocked atoms with the third level occupied. 



4-2. Particle-particle representation 

Another possible realization of the SU(3) algebra is: 



K21m 


-c f c f 


Ki2 TO 


^Q m C^ m 


Kn m 


= 1 — no m 


K3i m 




Ki3 m 


^Q m C^ m 


K22m 


= n lrn 


K32m 


= C \m C l m 


K23m 


= C lm C T"i 


K33 m 


= nf m , 



(48) 

where c£. m and c CTm are, respectively, fermionic creation and annihilation operators, 
and n CTm are the number operators cJ. m c CTm . The two lowering operators K21 and K31 
create a pair of particles and thus form a spinorial pair: 



K21m 



(49) 



The other lowering operator (K 32m = T +m ) is the ladder operator of a SUt{2) 
subalgebra of SU(3): 

[T +m ,P{j=Pj m . (50) 

It is well known that the particle-particle representation of the RG models associated 
with the SU(2) algebra describes the interaction between scalar pairs [T^], whereas the 
SO (5) version describes the vectorial pairing, i.e., the interaction between pairs forming 
a triplet [TS]. The SU(3) version presented here is associated with an interaction 
between pairs forming a doublet, we call it spinorial pairing or T = 1/2 pairing. As 
it was done in [15] for the SO(5) algebra, from a linear combination of the rational 
integrals of motion (|33p . the following pairing Hamiltonian for spinorial pairs results: 

H T =i/2 =^2 z m (n\m + n lm ) + g^2P} n ■ P n , (51) 



run 



where the scalar product is X} CT =T1 amF <*■"■> anc ^ we nave fixed the parameters 
^ 1 = 1/g and ^ 2 = 0. In this particle-particle representation of the SU(3) algebra, the 
Dynkin labels of the highest weight states are associated with the number of unpaired 
particles in each single particle level m and their transformation properties under the 
subgroup SUt(2). The number of variables Ek and tui is, respectively, the number of 
pairs and the number of |-particles in those pairs. 



4-3. Tavis-Cummings models for three level systems and two bosons 

In this section we derive another integrable model, which consists of L copies of the 
SU(3) algebra interacting with two different bosons. This model is derived from the 
trigonometric version of the previous section in two steps: first we take the bosonic 
mapping of one copy of the SU(3) algebra (let us say the copy 0), and then we let the 
degeneracy fio of this bosonized copy go to infinite. 
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We assume the SU(3) irrep (AJ,Aq) = (O ,0) for the copy to be bosonized. In 
this case the bosonic mapping for the SU(3) algebra is [25] : 

K„ -► n - a f a - b f b K 22 -> a f a K 33 -► b*b (52) 

K 3 i -► bVflo - ata- btb K 32 -> b f a K 2i -> a f v^o - at a - btb, 

where a^ and b^ are boson operators. This mapping is hermitic, therefore the mapping 
of the operators K Q ^ with /3 > a can be obtained by conjugating the previous ones. 
In the limit fin — > oo the previous mapping is reduced to: 

Kn^flo-ata-b'b K 22 -> a f a K 33 -> b*b 

Kax^iv/f^bt Kga^bta K 21 i v / f^a t , (53) 

where the imaginary factor i is introduced for future convenience. We consider now 
the trigonometric integrals of the RG model (|3"4"]) for a system of L + 1 copies of the 
Lie Algebra: m — 0, 1, L. Then we bosonize the copy with the "0" label using (f53|) . 
In order to avoid divergences in taking the limit Oo — * oo, we define a new set of 
variables e m by using the freedom we have to choose the parameters z m . These new 
variables are defined through: 

cot(z m - z ) = -7= m = l,...,L (54) 
V "0 

From this definition and the trigonometric identity (cot(A — B) = (cot (A) cot(-B) + 
l)/(cot(_B) — cot(A)), it is easy to show that in the limit fio —> 00: 

cot(z m - z n ) -> — . (55) 

By substituting these results in the trigonometric integrals Ro and H m (|34[) . we obtain: 
K = -C X a f a- (C 1 +C 2 )b f b + C l 

L 

+ (a f Ki 2n + aK 2 i„ + bK 3 i„ + ^Kisn + e n K Un ) (56) 

n m = -c 1 K 22m -(c 1 + c 2 )K 33m + c m i+ £ ^ K °fr" K »" 

- (a f Ki 2m + aK 2 i m + bK 3 i m + b^Ki 3m + e m Kn m ) , 

where we have rescaled the operators Ro and R m (JZq = Ro/\/fio and lZ m = 
Rm/V5)i dropped the terms with inverse powers of V^o, and defined new variables: 

The constants in the integrals of motion ([56]) are Co = —(1/3) J2i £ ( n i + (2C 1 +C 2 )(^o/3) 
and C m = (e m n m )/3 + J2'i n l/ ( e * - e «0 + (2C 1 + C 2 )?W3. 

To obtain the eigenvalues of the previous operators, we take the eigenvalues (|43|). 
with Z(zj — Zi) and Z{Ek — Zi) substituted by cotangent functions, and then the limit 
fio — > 00. Before performing this limit we introduce, as we did in (|54[) . two new sets 
of variables associated with the parameters E/~ and u>k'- 

cot{z Q - E k ) = -—7== cot(z -u){) = — -=!=. (58) 
vi in vfin 
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The same trigonometric identity that yields (|55[) allows us to write the limit f2 — > 00 
of the following functions: 

, / ci \ V^o , / n V^h , / \ \A^o 

COt(2 m - hk) — > COt(2 TO - CJ;J — > COt(u>;< - CJ;J 



- £m j/j - e m yi - yv 

cot(ff fe - -> ^ cot(-E fc - 03) -» — , (59) 

where we have used the definition ([54")) . With the limit values of the previous cotangent 
functions, the eigenvalues ([43]) become: 

L Mi 

T = ^2 k ln£n +Co-^2%k (60) 



k=l 



r m = -C 1 fc 2m -(C 1 +C 2 )fc 3m +C m + E3a=lKmkan ~ki m e m 

Mi 1 M 2 ^ 

+ (k lm -k 2m )^ Z K^m-fem)^] ~ , (61) 

k—1 l — l 

which are the eigenvalues of the operators ([55)1 . Note that the term Co, that diverges 
in the limit f^o — > oo, appears explicitly both in the operator and its eigenvalue, then 
we can easily get rid of it. e m are free parameters and the variables x k and yi are 
solutions of the Richardson-Bethe equations in the limit £lo —f oo. To derive them, 
we consider the Richardson-Bethe equations in the trigonometric case, which can be 
read from equations (|44p and ([4"5|) with Z substituted by cotangent functions. Then, 
using ([59]) we obtain: 

Mi q m 2 1 l , , 

k'(k'jtk) 1=1 m=l 

Ml 1 M 2 „ L , , 

-E^Er^-ET^-C 2 (/=1,.,M 2 ). (62) 

fc=l Xk yi V{ v^i) yi ' yi m=l £m ^ 

Note that the resulting Richardson-Bethe equations are identical to those in the 
rational case except by the linear term in the hrst line. The integrals of motion 
([35]) are now: 

L L 

N 3 = E K 33m + n b N 2 = E K 22m + n a 

m— 1 m— 1 

L 

Ni = E K n™ + °o - n a - n b . (63) 

m— 1 

The physical meaning of the Dynking labels (A^, A^) and n m is not modified by the 
introduction of the bosons, and it is the same already discussed above. The number 
of parameters x k and yi [M\ and Mq) in the Richardson-Bethe equations ([6"2"]) can be 
obtained from the general formula (I46|) . by extending the sum to the bosonized copy, 
and noting that the Young labels associated with it are [fcio, ^201 ^30] = [^Oj 0, 0]. 

One of the simplest Hamiltonian that can be derived from the previous integrals of 
motion, is obtained by considering just one non-bosonized copy of the SU(3) algebra. 
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Therefore L — 1 and we have 2 integrals of motion: R and R x . By taking a linear 
combination of these integrals we arrive to: 

H = z (Uo - C l) + z x (n x - Cil) 

= w a n a + uj b n b + 772K22 + 773K33 + g (a T K 12 + aK 2i + bK 3i + t^K^) , (64) 

where K Q ^ are the generators of the non-bosonized copy. The parameters of the 
Hamiltonian are related to the RG parameters by: 

x Lu a (A n - A w ) 2 A,, - A w u a \ - A w 7? 2 
C = T C = — ■ ei = -7 > 65 

g&u, g g&u 

and the constants zo and z\ are: 

z ° = ~g A A " A Zl = ~ g A A "a ' (66) 

with A^ = u)(, — uj a , and A v = 773 — 772 . The previous Hamiltonian describes the 
interaction between a set of Sli identical three level atoms and two modes of the 
electric held in the so called Rotating- Wave- Approximation [17] , and is a generalized 
version of the Tavis-Cummings model for three level atoms. 

If we write the Hamiltonian (|64|) in the particle-particle representation, we obtain 
the interaction between a doublet of pairs and a doublet of bosons: 

(P f ■ b + P ■ bt ) , (67) 

where we have redefined the bosons — > b_, bt — > b^, and the scalar product, -b, 
is P|b + + Pjb_. More complex Hamiltonians can be derived from the integrals 
but we will not discuss them here. 



5. Conclusions 



In this contribution we introduced a Gaudin algebra valid for any simple Lie algebra 
in the rational, trigonometric, and hyperbolic cases. With this algebra, we derived 
the complete set of quadratic RG integrals of motion, and found their respective 
eigenvalues. Focusing this formalism on the rank-two SU(3) algebra, we worked out 
in detail the RG models. For the rational case, expressions for the integrals of motion 
of degree three were obtained, and it was verified that their number is the necessary to 
match the number of integrals of motion with the number of degrees of freedom of the 
quantum model. Some physical applications of the SU(3) models were discussed. The 
physical meaning of the variables in the RG formalism were discussed for two different 
representations of the SU(3) algebra, namely, the particle-hole representation and the 
particle-particle one. By taking the trigonometric version of the models, we derived 
a new family of integrable models which are related to a generalization of the Tavis- 
Cummings model to three level atoms and two bosons. It was out of the scope of 
this contribution to explore all the possible physical applications of the SU(3) RG 
models, some examples were given and others are expected to come. We think that 
we have put the ground for more detailed studies, such as numerical studies of the 
solutions and comparison with approximative techniques in physical scenarios beyond 
the limits of traditional diagonalization methods. The formalism presented here can 
be easily applied to other Lie algebras. Some physically relevant examples are the 
SU(4) and SO(5) versions in the description of High temperature superconductors 
[24], the SO (8) version for isovector-isoscalar pairing [25], and the SU(n) versions 
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which include generalized Tavis-Cummings models for n-levels atoms interacting with 
n—l different bosonic modes. The formalism can likewise be extended to the elliptic 
RG models, where no linear term in the quadratic integrals of motion is allowed. 
Another interesting extension of the present contribution is to derive more general 
solutions to the Gaudin conditions (fT0|) . using the very well studied solutions to the 
Classical Yang-Baxter Equations [26]. Finally, the higher rank RG models can be 
useful to shed some light in the non completely well established definition of number 
of degrees of freedom in finite Quantum systems. In this contribution, the definition 
of intcgrability coming from the Yang-Baxter Equation was linked to the one coming 
from the equality between the number of quantum degrees of freedom and the number 
of independent integrals of motion. This latter definition, initially supposed to be 
exclusively related to dynamical symmetric models, can be extended, as shown in this 
contribution, to the RG models, which are not necessarily dynamical symmetric. 
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